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ABSTRACT 


The  various  components  of  the  COPSE  Acquisition  System  are  briefly 
summarized.  The  system  is  theoretically  analysed  to  determine  what  effect 
added  Gaussian  noise  will  have  on  its  performance.  On  linear  mode,  the 
Gaussian  noise  was  found  to  negate  the  quantization  effect  of  the  system 
provided  the  variance  of  the  noise  was  high  enough.  On  log  mode,  the  system 
was  found  to  perform  better  in  the  low  photo-diode  current  range. 


resume' 


Les  divers  elements  du  systeme  de  saisie  des  donnees  COPSE  font 
d'objet  d'une  breve  description.  Le  systeme  est  analyse  du  point  de  vue 
theorique  pour  determiner  les  effets  possibles  du  bruit  gaussien  ajoute  sur 
son  fonctionnement.  En  mode  lineaire,  on  a  constate  que  le  bruit  gaussien 
annulait  1'effet  de  quantification  du  systeme  lorsque  la  variance  du  bruit 
etait  suffisamment  elevee.  En  mode  logarithmique,  on  a  constate  que  le 
systeme  avait  un  meilleur  rendement  lorsque  la  photo-diode  avait  un  faible 
courant. 
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1.0  INTRODUCTION 


The  Coarse  Optical  Power  Spectrum  Estimator  (COPSE)  at  Defence 
Research  Establishment  Ottawa  consists  of  a  coherent  optical  Fourier 
transform  system,  a  64  element  solid  state  detector,  electronic  hardware  and 
software  run  on  a  NOVA  1200  and  a  XEROX  Sigma  9  computer  (Fig.  1).  The 
signals  coming  from  the  64  photo-diodes  in  the  detector  are  first 
multiplexed  onto  one  line.  The  multiplexed  signal  is  then  amplified,  biased 
and  amplified  again  before  it  is  sampled  and  digitized  by  an  A/D  converter. 
The  samples  are  subsequently  stored  on  magnetic  tape. 


The  various  components  of  the  COPSE  Data  Acquisition  System  were 
discussed  in  depth  in  an  earlier  report  [1].  Only  a  brief  summary  will  be 
presented  here.  Instead,  this  report  is  a  continuation  of  that  earlier  one. 
It  contains  a  theoretical  analysis  of  the  effect  that  additive  Gaussian 
noise  has  on  the  system. 


2.0  COPSE  DATA  ACQUISITION  SYSTEM 


The  data  acquisition  segment  of  the  COPSE  system  consists  of  a  64 
element  solid  state  detector,  electronic  hardware,  software  and  negnetic 
tape  storage.  The  electronic  hardware  is  further  subdivided  into  three 
separate  devices:  a  multiplex,  an  interface  and  an  8  bit  A/D  converter.  The 
signals  coming  from  the  64  photo-diodes  of  the  detector  are  first  switched 
onto  one  line  by  the  multiplex.  The  multiplex  allows  for  either  linear  or 
log  mode  of  operation.  Letting  i^  and  Vm(k)  represent  the  output  current  of 
the  kth  photo-diode  and  its  corresponding  multiplex  output  voltage 
respectively, 


vm(k)  -  -  |ik|  x  104  volts 


(1) 


on  linear  mode  provided 


-15  nanoamps  <  ifc  <  -1.3  milliamps 


and 


Vk>  ■  -10  -  2  '°SlO  I'k 


(2) 


on  log  mode  provided 


-3  nanoamps  <  i^  <  -1  milliamps. 


|i|J  is  the  absolute  value  of  ik  and  k,  the  photo-diode  element  number,  is 
an  Integer  between  1  and  64. 

The  interface  next  biases  and  amplifies  the  multiplexed  signal.  The 
device  gives  the  user  the  choice  of  several  gains  (G),  scaling  factors  (H) 
and  voltage  offsets  (Voff).  The  interface  output  voltage  takes  the  form  of 


Vj(k)  »  -  G[M  Vm(k)  +  Voff]  volts 


(3) 


4 


on  linear  mode  and 


V r (k )  =  -  G[Vm(k )  +  Voff  -8]  volts  (4 

on  log  mode  when  the  input  voltage  is  Vm(k ) - 

.  Substituting  equation  (1)  in  (3)  and  (2)  in  (4),  and  solving  for 
|ik]  in  terms  of  V j (k ) , 


.  VT  k)  ,  Voff 

iJ  =  —  x  10-“  + - 

kl  GM  M 


V, (k )  Vof 

Kl  ■  ioT(T  '  9  *T 


for  the  two  modes.  Hence,  an  incremental  change  in  |ik|  (a  | ik | )  will  result 
in  a  corresponding  change  in  V j (k )  of 

AVj(k)  =  GM  x  104  a  |ik|  (7) 


on  linear  and 


AVT(k)  = 


2G  4  Kl 

1"  1°  W 


on  log.  Linear  mode  gives  a  constant  resolution  throughout  the  photo-diode 
current  range  but,  on  the  log  mode,  resolution  depends  on  | i k  | .  The  two 

modes  have  the  same  resolution  in  situations  where 


|lk|  -lxJ0±  I  . 

1  ki  In  10  M 


Equations  (7)  and  (8)  are  plotted  in  Figure  2. 


fLm 
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The  A/D  converter  has  an  input  voltage  range  of  0  to  10  volts.  It 

converts  V  (k)  into  a  binary  number  whose  decimal  equivalent  U  satisfies 
I  * 


(Lk  -  0.5)S  <  V j (k )  <  (Lk  +  0. 5)S 


and  is  an  integer  between  0  and  255.  The  step  size  of  the  converter  in  the 
COPSE  system 


S  =  -1-  . 


If  the  system  is  noise-free,  the  COPSE  Data  Processing  program  will 
approximate  Vj (k )  by 


Vk>  =  Lks 


with  an  error  due  to  quantization  of 


eq  '  Vk>  -  »Iq  (k)- 

Note,  eq  varies  between  +0.5S  and  -0.5S. 

The  program  calculates  an  estimate  of  |ikJ  (|ik|  )  of 


K-  I'kI  -  x  10- 


on  linear  mode  and 


<kle  ■  kl  i(>eq/GH 


on  log  mode. 


3.0  NOISE  ANALYSIS 


To  include  noise  in  the  analysis,  the  COPSE  Data  Acquisition  System 
is  assumed  to  contain  a  series  of  strategically-placed  current  and  voltage 
zero-mean,  Gaussian  noise  generators.  Referring  to  Figure  3,  there  is  a 
current  generator  at  the  output  of  each  of  the  64  photo-diodes  of  the 
detector  and  voltage  generators  at  the  outputs  of  the  three  electronic 
devices.  The  noise  terms  are  further  assumed  to  be  additive  and 
statistically  independent  so  equations  (1)  through  (4)  now  take  the  form  of 


vj|(k)  - 

-|ik|  x  104  +  nk  x  104  +  Np|  volts 

(16) 

vjj(k)  = 

-  10  -  log10(nk  +  |ik|)  +1^  volts 

(17) 

Vj  (k)  = 

-  G[M  vjj(k)  +  Voff ]  +  Nj  volts 

(18) 

Yj (k )  = 

-  G  [vjj(k)  +  Voff  -  8]  +  Nj  volts 

(19) 

where 

nk*°k2  the  k**1  photo-diode  noise  current  and  its  variance, 

Nm* no^se  voltage  added  to  the  output  of  the  multiplex  and  its 
variance, 

Nj,oj2  noise  voltage  added  to  the  output  of  the  Interface  and  its 
variance. 


( 


KJ 

1 

& 


The  superscript  N  on  V^(k)  and  on  vjj*(k)  denotes  a  signal  which  has  been 
corrupted  with  noise. 


The 

equivalent 


decimal 


A/D  converter  converts  Vj(n)  into  a  binary  code  whose 
,  represented  now  by  L^,  must  satisfy 

(l{!  -  0. 5)S  <  vj{k)  <  (Lj|!  +  0. 5)S. 


Lj*  is  an  integer  between  0  and  255.  Hence,  the  error  due  to  quantization 


e“  -  <k>  -  Vjq (k) 


V?q(k) 


Before  the  digitized  signal  is  processed,  it  is  once  more  corrupted  with 
added  noise  at  the  output  of  the  interface. 


Lk  +  Lk  +  Ma> 


Na  being  the  noise  voltage. 

The  interface  output  for  each  photodiode  is  sampled  and  digitized  64 
times.  To  average  out  the  noise,  the  processing  program  calculates  the 
averages  of  these  samples  and  estimates  Ji^]  to  be: 


<kL  *  kl  -  "k  -  l"i  - ♦  vl 


on  linear  mode  and 


kl,  ■  kl  to 


on  log  mode. 


^  (Na  s  +  "i 


I*  "  J  ToQio t1  +  ITT] 

in  2  in  Pkl 


•  VJ 


•*>  -*•  ».  '  .  vW  V  *.**••*.*  V  V  V 


«S 


As  the  number  of  samples  taken  of  V j (k )  Increases,  the  averages  in 

(24)  and  (25)  approach  their  ensemble  or  expectation  values.  The 
expectation  value  of  a  Gaussian,  zero-mean,  random  variable  is  zero  [2]. 
Equation  24  therefore  reduces  to 


f  I  -  If  I.  iielLlor. 

">e  1  nl  GM 


The  notation  E{  }  stands  for  the  operation  of  ensemble  averaging.  If  all 
the  samples  of  Vj(k)  are  within  the  range  of  the  converter,  it  can  be  easily 

proved  that 


eq  +  0.5S  x  E{Vjq(k)}  -  VIq(k)  <  eq  -  0.5S 


where 


E(»^(k)|  -  Vjq(k)  =  eq  -  Efe"), 


E{vjq}  has  the  same  range  of  values  as  e_;  I.e.,  ±0.5S.  Futhermore,  the 
expectation  value  of  Vjq(k)  is  within  one  step  size  of  its  noise-free 

counterpart.  In  other  words, 

lE{VIq(k>}  -  VXq(k)  |  <  S.  (2$ 


A  more  rigorous  discussion  of  the  problem  is  contained  in  Korn's  book 
[3].  Using  his  results  and  after  some  mathematical  manipulation.  It  can  be 
shown  that  E{eq}  is  given  by  the  sum  of  weighted  sines: 

_  r  Mi  C  ft  Lilt  .  r*e„  2 IT,  ,  2  a  2 IT2*2  _ 


E{ejl-S  l  fsi«mx  exp  { 
q  n  i  =  1  4  s 


§2 - 1.  OO) 


Furthermore,  since  the  noise  terms  are  additive  and  statistically 
independent. 

Op2  a  oj2  +  G2  M2  x  10®  2  +  G2M20|y|2  (31) 

Obviously,  E{eq}  depends  on  the  step  sjze  of  the  converter,  eq  and 
the  noise  statistics.  The  larger  the  ratio  ^  is,  the  smaller  E{eN}  will  be 
and  the  better  the  estimate  of  |ik|  will  be.  Higher  gains,  scaling  factors 
and  noise  variances  will  accomplish  that.  See  Figure  5. 


S  <  3on, 


only  the  first  term  of  the  summation  in  (30)  contributes  appreciably  to  the 
estimate.  Equation  (30)  can  be  approximated  by 


,2w  e„,  -2*2  a. 


E(e  "I  *  1  sin  1^)  exp 


More  terms  must  be  included  In  the  summation  if  S  is  larger  or  on  is 
smaller.  Nevertheless,  even  when 


S 


(34) 


E{ejj}  is  less  than  its  corresponding  noise-free  value  (Table  1).  The  noise 

tends  to  negate  the  effect  of  quantization  provided  the  noise  variance  is 
high  enough.  In  fact,  in  situations  where  eQ  is  zero  or  ±Q.5S,  it  is 
theoretically  possible  to  recover  | 1 ^ |  exactry. 

Experiments  done  with  the  COPSE  system  indicate  that  Vjq(k)  has  a 

range  of  1  to  6  step  sizes  depending  on  the  element  number.  The  higher  ring 
elements,  due  to  their  larger  areas,  have  noiser  output  currents. 

It  is  a  common  practice  in  data  processing  schemes  to  correct  for 
dark  current.  In  the  COPSE  system,  the  correction  can  only  be  done  when  the 
system  is  on  linear  and  there  is  no  voltage  offset.  The  correction  was 
found  to  be  insignificant. 


3.2  Log  Mode  of  Operation 


Equation  (25)  becomes 


(kle  -  kl 


10 


-  E(eg} 
2G 


10 


E{log10(l  +  jy-j-)} 


(35) 


in  the  limit  as  the  number  of  samples  taken  of  Vj(k)  becomes  large. 
Assuming  that 


-1  < 


1, 


(36) 


which  is  a  reasonable  assumption  to  make, 


±  S/16 

(±  6.25  x  10"2S) 

±  S/8 

(±  1.25  x  lO^S) 

±  3S/16 

(±  1.875  x  lO^S) 
±  S/4 

(±  2.5  x  10-lS) 

±  5S/16 


S  =  2on 

S  =  3  on 

0 

0 

8.761  x  10"3S 

1.357  x  10-2S 

1.619  x  10-2S 

2.508  x  10"2S 

2.115  x  10"2S 

3.281  x  10-2S 

2.289  x  10- 2S 

3.551  x  10- 2S 

2.115  x  10'2S 

3.282  x  10“2S 

1.619  x  10-2S 

2.513  x  10"2S 

8.761  x  10"3S 

1.361  x  10-2S 

0 

0 

the  following  approximation  can  be  used: 


lo9  10 


(1  *  h] 


_i _ fnk 

In  10  ‘JlTT 


(37) 


and,  since  nk  is  a  Gaussian,  zero-mean,  random  variable,  equation  (35) 
becomes. 


*k 


10 


°k2  J _ 

|ik[2  In  10 


(38) 


|ik|  is  not  only  limited  by  E {e^}  but  also  by  the  statistics  of  the  photo¬ 
diode  noise  current  and  |ik|2.  The  lower  the  noise-free  photo-diode 
current,  the  better  will  be  the  estimate  of  |ik|.  See  figure  4. 


If  all  the  samples  of  V^(k)  are  within  the  range  of  the  A/D  converter. 


-0.5S 


G  ok2 

+  r"- 10  KF 


v j (k)  -  E  {V]fq (k )}  <  0.5S  + 


—  y2 

In  10  • 


(39) 


Consequently,  E(Vjq(k)}  is  within  one  step  size  of  Vjq(k)  even  though  its 
range  of  possible  values  are  biased,  i.e.. 


-S  + 


in  10  jtj*  '  Wk)  '  E,vlq(k)l  ‘  S  *  T^io  frk7 


(40) 


Again  referring  to  Korn's  book  and  after  a  great  deal  of  mathematical 
manipulation.  It  was  found  that 

E*ei*  " » *  1 1  (i  +\2k2 )'’  expl"fr  *2  °n«2!  x 

Sin  (if  e,  -  t.n-1  (K«)  *  f  L  ^,1  (41) 
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where 


4Gn  o^2 

K  ■  sVTo  Kfa 


°ru2  = 


?  9  °k  1 

G  *  °i  *  wTo  Kp-  rrw  • 


°ne2  1S  Plotted  In  Figure  6  as  a  function  of  ry and  i.  Note,  o2.  exhibits 


an  exponential-like  decty  with  i  and  is  a  maximum  at 


KF 


S  In  10 


Furthermore,  since  each  term  in  the  summation  in  (41)  is  multi  pled  by  a 
function  that  goes  as  the  inverse  of  A2,  one  would  not  exoect  the  higher 
terms  in  the  summation  to  contribute  significantly  to  E {e^} . 

In  the  case  of  very  weak  photo-diode  noise: 


V 

KF 


«  i , 


equation  (41)  reduces  to  a  form  similar  to  that  of  equation  33: 

r  m.  S  °?  f2nA  .  -2irA2 

Eleq)  *  -  t  l  !  sfn  I—  eq)  exP  (— j—  “n2) 


and  (see  Figure  6) 


an2  =  G2  au2  +  or2  +  — - - 

"  M  1  In2  10 


weak  noise;  slope  =  4G2/ln210 


C,  =  16  iri  ;  C2  =  256r»l< 

S  In2  10  S«ln«  10 


I 


rj 


4Girl 
S  In  10 


a 2  +  4G2  <r 


M  I 


In2 10 


£k  _ 1 

W  1  +  1 


;  K  *  4Gir  at 
Sin  10 


Fig.  6:  A)  Plot  of  on»2  Versus  tttt 

Pkl 

B)  Plot  of  on,2  Versus  i  on  Log  Mode 


(S  *  step  size  of  converter;  G  =  gain;  oM2,  or2,  ok2  *  noise  variances) 
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If  e_  is  zero  or  ±0.5S,  Efeq}  is  zero  but,  according  to  equation  38,  | ik  j e 
is  not  equal  to  |  i  k  | .  It  is 


r*k~F  10 


> 


(48) 


The  biasing  effect  of  the  log  mode  does  not  allow  for  this.  Complete 
recovering  of  |ik|  occurs  only  if 


E 


-6  V 

In  10  |ik|2 


(49) 


Experiments  done  with  the  COPSE  system  on  log  mode  indicated  that  j 

Vjq (k )  has  a  range  of  values  of  1  to  3  step  sizes.  The  lower  photo-diode  ^ 

currents  have  a  wider  range  of  sample  values  as  expected. 

k 

k 

V 
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